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Abstract 

In this paper we discuss the presence of temperature-dependent squeezing in the collective exci- 
tations of trapped Bose-Einstein condensates, based on a recent theory of quasiparticle damping. 
A new scheme to measure temperature below the critical temperature (T <C T c ) is also considered. 
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I. INTRODUCTION 



Following the first experimental observation of Bose-Einstein condensation in 1995 1], 
analogies between this phenomenon and non-linear optics have been pointed out (3, 3|. 



Recently, squeezing of the matter waves has been analysed |4|, |5|, |6[; by considering the 
role of the nonlinear interactions between the atoms in generating entangled atomic beams, 
namely the spin-exchanging collisions of spinor Bose-Einstein condensates (BEC's). 

Here we analyse squeezing in relation to the coupling of collective excitations in a Bose 
condensed gas in a trap. Most of the observed effects of elementary excitations can be ex- 
plained using the standard Bogoliubov theory However, it cannot account for higher 
order processes such as the Beliaev damping of the excitations observed recently for a scis- 
sors mode of a BEC in 87 Rb 0. Beliaev damping is a mechanism in which a quasiparticle 
with energy E 2 interacts with the condensate and generates two quasiparticles with energy 
Ei. The process is actually based on a four- wave process that is desguised as a three- wave 
interaction analogous to degenerate parametric down-conversion. The new quasiparticles, 
being produced in pairs, constitute a squeezed twin phonon beam. Recent theoretical de- 
velopment including the presence of interactions between the quasiparticles f9| provides an 
explanation for such processes. Finding quantitative estimates for the amount of squeezing 
and entanglement based on this theory is one of the main tasks of this paper. 

Considering a system containing 100 87 Rb atoms in a spherical trap, we calculate 
the degree of squeezing of the quasiparticle excitations. A related experimental challenge 
that we also address is the accurate determination of temperature in BEC's at T <C T c . 
The temperature is currently estimated by fitting a thermal Gaussian profile to the atomic 
cloud. We propose a different way to determine temperature by observing the variation of 
the envelope of collective oscillations. 

The paper is structured as follows: In Section |Hj the Hamiltonian in both particle and 
quasiparticle basis is discussed in relation to squeezing. The Hamiltonian in the quasiparticle 
basis accommodates processes such as the Landau and Beliaev damping. In Section ILTT| we 
obtain equations of motion for the quasiparticle modes and give quantitative estimates for 
the degree of squeezing and entanglement in the quasiparticle excitations. The expected 
"damping rate" at various temperatures is calculated. This is then used as a temperature 
measurement calibration for BEC. In Section Hv| we conclude. 
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II. MANY-BODY HAMILTON! AN 



Within the second quantised formalism, the many-body Hamiltonian for a system of 
bosons with binary interactions can be written as, 

H = H-?a\aj + - ^(ij\V\km)a\aja k a m , (1) 

ij ijkm 

where the matrix elements H*? are given by 

Hg = f d 3 r^(r)^i(r). (2) 

Here, H sp = — <^V 2 + V trap , is the single-particle Hamiltonian with a confining potential 
Vt r a P , and the basis state wave functions are ipi(v). The matrix element (ij\V\km) denotes 
the matrix element for a bare interaction potential V(r) between two atoms; in this work, 
we shall consider the typical delta-function contact interaction potential for V"(r). The 
operators a\ and dj are the creation and annihilation operators for mode i that obey the 
Boson commutation relations 

a]] = 5ij, [Sj, %•] = [of, at] = 0. 

The Hamiltonian (0) is written in a single-particle basis where the operator dj annihilates 
a particle from the state with wave function ipi(r). The wave function ^o( r ) describes the 
condensate, whilst the remaining functions form a complete set orthogonal to the condensate. 

A. Particle Hamiltonian 

Let us first investigate the condensate by concentrating only on the condensate mode. 
This one-mode approximation applies well for temperatures T ~ which is indeed the case 
for many experiments. We find the Hamiltonian is then simplified to 

H = e 00 ala + ^(00|V|00) aJaJa ao- (3) 

Treating the deviations from the ideal condensate as perturbations, we are able to expand 
a as 

d = z + c, (4) 

where z = ^/Nq represents the mean-field (condensate), and c is the small quantum per- 
turbation, c may take the meaning of non-condensed atoms. We choose the phase so that 
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z is real. The mean-field approximation, although not exact, has been very successful in 
describing a large number of experimental results. Using the expansion for So and splitting 
our Hamiltonian into a Gaussian and a non-Gaussian parts, we get the Gaussian part of our 
Hamiltonian as 

H = e [z 2 + z{c + gt) + g f g] + (00| ^ |00) [z 4 + z 2 (c 2 + c t2 ) + 2z 3 {c + gt) + 4z 2 gtg] . (5) 

The non-Gaussian part may be treated using perturbation theory [ljj]. The resulting 
Heisenberg equations of motion for c and & are 

(6) 

where a = ± (e Q + 2(00\V\00)z 2 ), b = ±(00| V\00)z 2 , and k = J (e z + (00\V\00)z 3 ). The 
solutions of these equations are 

i A r \ / z(r\\ \ / n I 

(7) 

where A = cos Tt — i 4 sin Tt, i? = — i ^ sin Tt, D = cos Tt — i £ sin Tt + ^ b ^ K , and 







T = Va 2 — b 2 . It is important to note that the above equations are formally equivalent 
to dynamic equations that describe parametric down conversion in quantum optics ^j] if 
k = 0. 

Defining the quadrature operators as X p {t) = c ( t )+ ct ( t ) ari d Y p (t) = — ~ ct — , we have that 
the variances can be written as 

(AXMf = <^±I2 ^ Tt + ^! 8 m 2 Tt) (8) 

{ AY M y = <^ (cos* Tt + ^ sin* Tt j (9) 

where n is the number of excited particles. In this case, the magnitude of a and b, or 
equivalently the magnitude and phase of the mean field, z, control the amount of squeezing 
of the particule quadrature operators. 



B. Quasiparticle Hamiltonian 



On the other hand, it is possible to describe the system including the higher order pro- 
cesses. The dominant ones correspond to the quadratic Hamiltonian Hq — H q + Hi + H 2 , 
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where the subscripts 0,1,2 indicate the order in a*. This may be diagonalised exactly using 
the standard Bogoliubov transformations, changing the Hamiltonian to a quasiparticle basis 
where the quasiparticle operators $ are defined by 1 

Here and V^- are the well-known matrices associated with the Bogoliubov trans- 
formations, obeying the orthogonality and symmetry conditions UU^ — VV^ = 1 and 
UV T — VU T = 0. These matrices have to be evaluated numerically depending on tem- 
perature and the geometry of the trap. The non-quadratic terms are expected to be small 
and can be dealt with perturbatively. The energy and shape of the condensate change when 
one includes these terms. It has been shown |9( that in terms of the quasiparticle operator 
Pi the effective Hamiltonian can be written as 



H' = const+ S ^{e i +/\e i )p\p i +\ 



E [^kPiPjPk+djkPjPjPkl +J2viPi 



+ h.c. > , (11) 



where the const term simply shifts the zero of energy and Ae^ is the energy shift from first 
order perturbation theory. The Hamiltonian given by equation (jll[) contains terms beyond 
the Bogoliubov approximation, showing that we are treating interacting quasiparticles. We 
are thus taking into account important processes such as Landau and Beliaev damping. 
Landau processes, in which two quasiparticles collide to form a single quasiparticle, can 
not occur at zero temperature because there are no excited quasiparticles. However, they 
are dominant at high temperatures. Expressions for the coefficients of equation are 
given in Appendix El As we are interested in finite temperature, the U and V matrices 
are calculated self-consistently using the Bogoliubov-de Gennes (BdG) formalism for each 
temperature. From the Hamiltonian (fTTj) . the Heisenberg equation of motion for P p is 

P = \ + ^pPp + Yl + E Pi*PH + E v i*PkPi> ( l2 ) 

'" j,M0 j,k^0 



1 For simplicity from now on we drop the explicit notation in the operators 
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where 

Up = {e p + Aep)/h, 

"jk Cpjk/^y 

Pjk = {^pjk + £jhp + £kpj) l^i 

Vjk = (Cj pk + Cjkp) /h- 

Equation (|12j) and its complex conjugate are valid for the general case. However, a 
dominant process can be defined by choosing the appropriate geometry for the trap. In this 



form, the number of states involved is reduced [12|. The trap geometry can be modified by 
adjusting the frequencies in the radial and axial directions independently and the selection 
of a dominant mode in this way has been demonstrated jlflj . In particular, the observation 
of a Beliaev process has recently been reported for a scissors mode, where one mode is 
resonantly coupled to two modes of half the original frequency this is the kind of process 
that we seek to model in the present paper. We note that terms in £ and rj describe the 
spontaneous decay of quasiparticles and thus do not generally conserve energy. Therefore, 
we will only keep terms containing a and v. More details of these arguments can be found 
in reference and the main points are summarised in Appendix El For a Beliaev process, 
the equations of motion for modes p = 1 and p = 2 are 



Pi = -i^-m^l (13) 

ft = -io^-i^iftft. (14) 

(15) 

plus their complex conjugates. The factor of 2 that appears in equation (fTf)) comes from 
the fact that a\\ = v 2 i/2. 

Considering the system at a temperature T = 20 nK, with iV = 10000 particles and a 
trap frequency uit ra p = 2ir x 100 Hz with spherical geometry and using 87 Rb atoms for which 
the s-wave scattering length a = 10a , where a is the Bohr radius, the coefficient u 2 i equals 
1.90 x 10~ 2 in trap units. The equations of motion for {3% and (3\ are profoundly affected by 
temperature through the values of v 2 i- 

Equations of this sort have extensively been studied in quantum optics. In this case we 
apply them in order to study squeezing in the quasiparticle excitations and its temperature 
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dependence. It is important to mention that the damping of excitations has earlier been 
described in terms of nonlinear mixing this work differs from that description as the 
operator nature of the quasiparticle annihilation operator fa is retained, and the calculation 
is not restricted to the quadratic approximation. We are in effect extending the work of Ref. 
(if] ], such that spontaneous quantum processes are included. 

III. TEMPERATURE DEPENDENT COUPLING PROCESS 
A. Non-depleted regime 

In the case in which mode 2 has a much larger population than mode 1, we may approx- 
imate the operator 02 by a c- number 62 and ignore the depletion of mode 2. The solutions 
of the Heisenberg equations of motion for the operators fa and j3\ are then given by 



(16) 



fa(t) \ I cosh(fit) -isinh(ftt) \ I /3 X (0) 
p\(t) J \ isinh(m) cosh(fit) J \ ^(0) 
where Q = vi\ &2, an d we have chosen the phases such that the coefficients are real. 
This represents the physical situation where the pump is continuously driven by a resonant 
excitation. It is important to mention that the approximation is no longer valid when a 
noticeable down- conversion has occurred. 

Solutions (|16|) are equivalent to equations obtained in the description of degenerate para- 
metric down-conversion in quantum optics. This process is well-known to be an efficient 
source of squeezed states. In this case, the corresponding squeezing parameter is 



r 



fit, (17) 



which is temperature dependent via fi. 

Analogously, one can define the quadrature operators for (3 and ft and calculate the 
variances 



, 2 (2JVi + 1) 



(AX(t)Y = v L - y exp(-2fit), (18) 



(AF(t)) 2 = ( 2iV i + h exp ( 2 Qt) ; (19) 
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where N± is the number of particles in mode 1. In optics the quadrature operators are 
well-defined quantities corresponding to the amplitude and the phase of the electromagnetic 
(EM) oscillation. In a similar fashion, we give our quadrature operators X and Y the 
interpretations of amplitude and phase of oscillations. The dependence of equations (|18|) 
and (|19|) implies that the amount of squeezing is related to the number of lower-mode 
atoms present. Consequently, it is connected to the temperature of the ultracold atoms 
in the condensate defined in terms of the initial Bose-Einstein (BE) distribution of 
quasiparticles. 

However, in an experiment what one measures some correlation functions, for instance 
(f3\f3i). Assuming an initial number state in the quasiparticle basis, some important corre- 
lation functions are given by: 



A plot of the correlation function (|2T)|) is shown in Figure El a). A non-zero value for this 
correlation function implies the presence of squeezing and hence of entanglement between 
the two modes. The quantity given by equation (}2*Tj) carries information about the behaviour 
of the population in mode 1 and it is directly measured as the amplitude of oscillation. A 
plot of this correlation function as a function of time (t) and temperature (T) is given in 
Figure El b). At t = the behaviour is determined by the BE distribution of Ni. As time 
evolves the population of the lower mode increases. 

B. Depleted regime 

Let us now include the effect that depletion of mode 2 has in the description of the 
coupling process. In order to tackle the problem, we will consider the equations of motion for 
the number operators N{ = /3|/3j. They can be calculated through the Heisenberg equations 
1I1HI14I) . namely 




(20) 



(21) 




(22) 



(23) 
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It is possible to uncouple the equations (|22j) and (|23|) using the constant of motion 
A = Ni + 2N 2 and calculating the second derivative. The uncoupled equations are given by 



v 2 ! 1 (-3Nf + 2AN 1 + A), (24) 



d 2 Nx 
dt 2 

A 2 AT 7/ 2 

^ = ~f{l2Nl - 8AiV 2 + A" - A). (25) 

When the population of mode 2 is large, it is possible to determine the time evolution by 
taking the average of equation (J25j) . In that case we have a c-number second-order differential 
equation with the following initial conditions 

diV 2 



dt 



= 0, (26) 
N 2 (0) = N 2Q ; (27) 



The solution can be expressed as follows 



' Nw + ( aa _ iv 20 ) sn 2 (m^a^ ^g=jgj) , for iV 20 < a 2 , 
N 2 (t) = I (28) 

k «i + (iV 20 - ax)nd 2 (f^fi, ^/^f) , foriV 20 >a 2 , 
where sn(w, /c) and nd(ii, /c) are the standard Jacobi elliptic functions [18(, and the coefficient 
«i and a 2 are the roots of the quadratic polynomial P(N 2 ) = —N% + (A — N 20 )N 2 — N% — 
(A 2 — A) /A + N 2 qA. It has been already mentioned that the coefficient z/ 2 i is temperature 
dependent; therefore the population will change accordingly. The initial population was 
calculated for different temperatures as prescribed by the BE distribution. An initial driving 
was also taken into account. The solution was then calculated using these parameters; a 
plot is shown in Figure El as a function of temperature and time. 

The connection between the theory and experiment can be accomplished by observing 
that the initial part of the population surface can be fitted by a curve of the form A\ cos(rt). 
The parameter Y is related to the solution (|2*5|) as 



(29) 

l; 127V 20 1 ; 

Following this line of analysis, a fitting routine was used to obtain a plot of V as a function 
of temperature (FigureHJ). We can see from the plot that for lower temperatures the influence 
of the initial population is preeminent. For higher temperatures the coupling coefficient z/ 21 



dwindles quickly, thus the effect of the initial population is hidden. It can clearly be seen 
that there is a competition between the initial population and the coupling coefficient. The 
quantity T is notable in that it is a readily measurable temperature dependent quantity that 
links to the squeezing of quasiparticles. 

IV. DISCUSSION 

The detection of squeezed states in the case of the electromagnetic field can be achieved 
through well-known techniques, e.g. homodyne detection. On the other hand, scattering 
experiments have also been proposed as a way to detect squeezing in bosons [l^. The 
proposal makes use of a particle scattering off a boson field. Provided the particle is prepared 
in the proper input state, it can either absorb or emit a boson, and these two scattering 
processes can coherently interfere when the boson field is in a squeezed state. The change in 
the rate of scattering into a particular state gives a signature that the field has been squeezed. 
A schematic setup for the detection of squeezing is shown in Figure El An incoming beam 
of frequency uo_ is divided into two by a beam splitter, with one of the beams frequency 
shifted to u+, which may be achieved through acusto-optic modulation. The two beams 
are made incident into the boson field, e.g. BEC. Particles scattered off the squeezed field 
are detected by the counter. This would provide an experimentally feasible way to directly 
measure squeezing in BEC. 

In conclusion, we have described a process analogous to parametric down-conversion 
in trapped Bose-Einstein condensates. A fully second quantised order theory has been 
used to depict the condensate at finite temperature. The process corresponds to Beliaev 
damping, where a quantum of higher energy is divided into two quanta of lower energy. The 
coupling process is temperature dependent; this implies that any amount of squeezing in the 
elementary excitations of the condensate would give a direct indication of the temperature. 
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APPENDIX A 



The coefficients that appear in equation (|TT|) written in the position representation (using 
the contact potential) are determined by 



Cijk = VNqUo J d 3 r[^(r)w i (r)M j (r)u fc (r)+^o(r)^(r)^(r)t; fc (r)], (Al) 

Cijk = v^Vo^o J d 3 r{^*(r) [u*{v) U] {r)u k {r) + v*{v) Vj {r)u k {r) + v*{r) Uj {r)v k {r)) 

+ <Ao [< (r)M i (r)u fc (r) + u*(r)v j(r) u k (r) + v*(r)vj(r)v k (r)] } , (A2) 
Vi = ~ (Cm* + C qiq ) N q . (A3) 

where the indices i, j and k are labels that denote eigen-energy levels. For a 3-D condensate, 
the level i is an implicit notation for the quantum numbers n, I and m. 



APPENDIX B 



To prove that Beliaev damping is indeed the dominant process with the given trap ge- 
ometry discussed in the paper, it is only necessary to calculate various contributions from 
the framework of second order perturbation theory. The energy shift due to various higher 
order terms is given by the expression 

AE( P ) = -j2 9 , 'r*' J 1 + N *+ N A ( Bla ) 

+ ^jk±gi [1+JVi+Jv .] (Blb) 

+ J2 ^ ijP _ +Cipj ^ [N j -N i ), (Blc) 
ij^o '>■ ~ " + e i 

where line (jBla|) indicates spontaneous decay of quasiparticles. Here the coefficient is 
defined as a sum over permutations of the three indices in £ p j k , i.e. Cp^k = Cpjk + Cpkj + Cjpk + 
Cjkp + ikpj + Ckjp- Lines (jBlbjl and (jBlcj) give Beliaev and Landau damping processes respec- 
tively. Typical values are given in Ref. pjj and owing to the resonance in the denominator, 
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(|Blb|) may be made dominant with a correct choice of experimental parameters [1 
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Figure Captions 



1. Schematic diagram of energy levels for a BEC in a trap. In the Beliaev process, a 
quasiparticle of frequency 0J2 interacts with the ground state generating two quasipar- 
ticles of frequency uj\ that divide the initial energy equally. This process is analogous 
to the optical parametric down- conversion. 

2. a) Surface plot of the average A non-zero value of this quantity implies the 
presence of squeezing and hence entanglement, b) The average (PiPi) describes the 
behaviour of the population in mode 1. This solution is valid within the non-depleted 
regime. As discussed before, the behaviour at t = is determined by the BE distri- 
bution function. 

3. The population N2 as a function of temperature and time. 

4. Plot of T against temperature. This parameter is obtained by fitting a curve of the 
form Ai cos(Ti) corresponds to the initial part of the envelope of the population N 2 
shown in figure E3 

5. Schematic figure for the detection of squeezing by particle scattering. 



14 



2 



1 

Ground 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



■ T/T c 



Neutron 
source 




s 



